SUB-BARRIER FUSION OF DRIP-LINE NUCLEI 



K. HAGINO 
'^--v ' Institute for Nuclear Theory, Department of Physics, 

Q I University of Washington, Box 351550, Seattle, WA 98195, USA 

f~^ . E-mail: hagino@phys.washington.edu 

(N 

-._ A. VITTURI 

A\ • Dipartimento di Fisica, Universita di Padova and INFN, Padova, Ita 

r^ ' E-mail: andrea.vitturi@pd.infn.it 

^ ' • We discuss the role of break-up process of a loosely-bound projectile in subbarrier 

fusion reactions. Coupled-channels calculations are carried out for ^^Be -|- ^"^Pb 
and *'^He -|- 238^ reactions by discretizing in energy the particle continuum states. 
p^ , Our calculations show that the coupling to the break-up channel has two effects, 

C^ , namely the loss of flux and the dynamical modulation of fusion potential. Their 

net effects differ depending on the energy region. At energies above the Coulomb 
C 3 ' barrier, the former effect dominates over the latter and cross sections for complete 

^y\ ' fusion are hindered compared with the no coupling case. On the other hand, at 

CD ^ below the barrier, the latter effect is much larger than the former and complete 

f**^ , fusion cross sections are enhanced consequently. 
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1 Introduction 



^ ' Many questions concerning the effects of break-up process on subbarri£Lij_&j^ii£j] 

^ . have been raised-duxing the last few years both from the experimentall^'aHu'EH 

C ' and theoreticalQu'BllJ points of view. The issue has become especially relevant 

in recent years due to the increasing availability of radioactive beams. These 
often involve weakly-bound systems close to the drip lines for which the pos- 
sibility of projectile dissociation prior to or at the point of contact cannot be 
C^ ' ignored. 

One interesting question is whether the break-up process hinders or en- 
hances fusion cross sections. In addressing this, it is vital to specify which 
quantity is being compared with which quantity; otherwise a comparison is 
meaningless. The latter quantity is rather obvious in theoretical calculations, 
because one can artificially switch on and off couplings to the break-up chan- 
nel. When we refer in this paper to that break-up enhances/hinders fusion 
cross sections, it will be in this sense unless we mention explicitly. As for 
the former quantity, there are two possibilities, namely complete fusion cross 
sections or the sum of complete and incomplete fusion cross sections. From 
studies of fusion of stable nuclei where break-up process is not important, we 
have learned that any coupling of the relative motion of the colliding nuclei 
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to nuclear intrinsic excitations causes large enhancements of the fusion cross 
section at subbarrier energies over the predictions of a simple barrier penetra- 
tion model. lij It may not be difficult to imagine that the same thing happens 
to the break-up channel as well; cross sections for inclusive processes, i.e., the 
sum of complete and incomplete fusion cross sections would be enhanced by 
couplings to the break-up channel. On the other hand, one could also argue 
intuitively that break-up process removes a part of flux and thus cross sections 
for complete fusion would be hindered. As will be shown below, fusion cross 
sections are determined by the competition of these two mechanisms. 

In passing, this sort of consideration is relevant also when one discusses 
experimental data and compares them with theoretical calculations. For i: 
stance, it is important to beaiuin mind that the recent Padova/RIKEN data 
as well as the Canberra data,Q which used ^Be beams, contain both the com- 
plete and incomplete fusion cross sections for the neutron break-up channel 
(^Be — > n -f- *Be), while they are only the complete fusion cross sections with 
respect to the a- break-up channel (^Be -^ a + ^He or 2a + n). 

In the past, different theoretical approaches to the problem have led to 
controversial results, not only quantitatively but also qualitatively. Hussein et 
al. derived a local dynamical polarization potential Fdpp for break-up process 
and computed complete fusion cross sections asQ 

'^CF = ^E(2^ + l)4°^expQy'"dtW^(r(i))) , (1) 

where Pp is the fusion probability in the absence of break-up, k is the wave 
number in the entrance channel, and W{r) is the imaginary part of the dy- 
namical poralization potential T^dpp ■ The integral in the exponent is evaluated 
along the classical trajectory r{t), where r{t — 0) corresponds to the distance 
of closest approach. The idea of this formalism is that the effects of break-up 
can be well described in terms of the survival probability Psurv = 1 — Pbu = 
exp{2 J W{r) dt/h). Since W{r) is negative, the break-up process always hin- 
ders complete fusion cross sections in this formalism. Subsequently, Takigawa 
et al. pointed out that the classical trajectory used by Hussein et al. corre- 
sponded to the one for scatterinffii, and for fusion reactions the time integral 
should have been from — oo to O.El Consequently the effects of break-up were 
moderated over the estimate made by Hussein et al., but the qualitative con- 
clusion remained the same. These conclusions were later criticized by Dasso 
and Vitturi. They performed coupled-channels calculations by treating the 
break-up channels as a single state and obtained totally different results, i.e., 
enhancement of complete fusion cross sections due to the break-up. 13 



What would be the origin of these apparently controversial results? We 
argue that the coupling to the break-up channel leads to the dynamical modi- 
fication of fusion potential as in fusion of stable nuclei. At subbarrier energies, 
this effect is most relevant, leading to the enhancement of complete fusion 
cross section. Such effects are automatically included in the coupled-channels 
formalism. The dynamical modulation of fusion potential is related to the real 
part of a dynamical polarization potential, which both Hussein et al. and 
Takigawa et al. completely threw away. Therefore, if the real part of the 
polarization potential is properly taken into account, their formalism would 
provide an enhancement of complete fusion cross sections at energies below 
the Coulomb barrier. Actually, Ref. 12 shows that it is indeed the case, a' 
though their calculations are not satisfactory in a sense of Takigawa et c 
and a full calculation within their formalism has still been awaited. 

Another origin of the controversy might be the simplified assumption used 
by Dasso and Vitturi. As mentioned above, ihe entire continuum space was 
mocked up by a single discrete configuration, El and therefore the main feature 
of continuum couplings were not entirely included. This would underestimate 
the break-up effects, especially those effects, if any, which remove the incident 
fiux. 

From this point of view, in this paper, we repeat similar calculations as 
of Dasso arid, Vitturi, but by considerably increasing the number of contin- 
uum states. Ii3 We also use microscopic form factors in contrast to the previous 
studies, where form factors were assumed to have an exponential form. These 
two aspects enable us to describe simultaneously the dynamics of continuum 
couplings inside and outside the Coulomb barrier, leading to a more conclusive 
result on the effects of break-up process on subbarrier fusion. In the next sec- 
tion, we review briefly the coupled-channels formalism and define the complete 
and incomplete fusion cross sections. The specific applications of the formalism 
to '^"'^Be + ^°*Pb and ^'^He -I- ^^^U are given in the following section 3. 

2 Coupled-Channels Formalism for Subbarrier Fusion of Drip-Line 
Nuclei 

Our aim is to discuss effects of the break-up of the projectile nucleus on sub- 
barrier fusion reactions by solving coupled-channels equations. The-coupled- 
channels equations in the isocentrifugal approximation are given byll3 
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where the angular momentum of the relativ e rr iotinn in each channel has been 
replaced by the total angular momentum LI13 In eq. (||), V^ ' is the nuclear 
potential in the entrance channel, and e^ is the excitation energy of the n-th 
channel. Here, we assume that n = labels the ground state of the projectile 
nucleus and the all other n refer to particle continuum states, which are associ- 
ated with the break-up channels. It is straightforward to include bound excited 
states of the projectile and/or the excitations in the target nucleus. Fnm are 
the coupling form factors, which we compute on the microscopic basis. It is 
to fold the external nuclear and Coulomb fields with the proper single-particle 
transition densities, obtained by promoting the last weakly-bound nucleon to 
the continuum states. More explicitely, the single-particle form factor for the 
promotion from the bound state (ni^ijiTOi) to the continuum state {£2J2m2) 
with continuum energy E assumes the form 

= ^ ^(-1)™^+^ ,5(^1 -K £2 + A, even) /jilj2 - ^|A0 
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The functions REi2h{''') and Rnii-iji{r) are the single particle wave functions 
for the continuum and bound states, respectively. The potential Vr is the mean 
field felt by the single particle due to the presence of the target, responsible 
of the transition, not to be confused therefore with the projectile mean field 
generating the single particle wave functions. This potential obviuosly involves 
both a nuclear and a Coulomb component. Note that for long ranged transition 
densities, as in the case of a weakly bound system, the resulting Coulomb 
formfactor will differ from the pure r^^^^ form much outwards than in the 
traditional case of stable systems. See Ref. [16] for more details. 

The coupled-channels equations (0) are solved by imposing the incoming 
wave boundary condition (IWBC), where there are only .incoming waves at 
T^min which is taken somewhere inside the Coulomb barrier. EiltZI The boundary 
conditions are thus expressed as 

ipn{r) — > T„exp f -i / kn(r)dr' j r < r„iin, (4) 

-^ HJ;^\knr)Sn,0 + i?„i?/ + ^(fc„r) r > Vmax, (5) 



where 
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is the local wave number for the n-th channel and A:„ = \/2^{E — e„)/?i . Hi 

and Hi are the incoming and the outgoing Coulomb functions, respectively. 
The fusion probability is defined as the ratio between the flux inside the 
Coulomb barrier and the incident flux. For our boundary conditions given by 
eqs. (||) and (||), it reads 

fco 

for the n-th channel. Complete fusion is a process where all the nucleons of the 
projectile are captured by the target nucleus. We thus define cross sections pf 
complete fusion using the flux for the non-continuum channel (i.e., n—Q) asliS 

Kq ^ Kg ^ KO 

The flux for the particle continuum channel {n ^ 0) are associated with in- 
complete fusion, whose cross sections we define as 

3 Results 

S.l The "5e + ^osp^ reaction 

We now apply the coupled-channels formalism presented in the previous sec- 
tion to fusion of drip-line nuclei. We first consider the fusion reaction ^^Be + 
^"^Pb, where the projectile is generally regarded as a good example of a sin- 
gle neutron "halo". In a pure single-particle picture, the last neutron in ^^Be 
occupies the 2si/2 state, bound by about 500 keV. The strong concentration 
of strength at the continuum threshold evidenced in break-up reactionslij has 
been mainly ascribed to the promotion of this last bound neutrniL_tp the con- 
tinuum of P3/2 and pi/2 states at energy E^ via the dipole field. E3l3 Since the 
presence of the first excited lpi/2 state (still bound by about 180 keV) maji 
perturb the transition to the corresponding pi/2 states in the continuum, Ell 



we prefer here to consider only the contribution to the P3 ^ states. The cou- 
phngs to the pi /2 states are certainly not negligible, but we expect that they 
will simply further enhance the effects on the fusion cross section and therefore 
will not alter our qualitative conclusions. The initial 2si/2 bound state and the 
continuum P3 /2 states are generated by Woods-Saxon single-particle potentials 
whose depths have been adjusted to reproduce the correct binding energies for 
the lp3/2 and 2si/2 bound states. In particular, one needs for the latter case a 
potential which is much deeper than the "standard" one. As we mentioned in 
the previous section, the form factor F{r; Ec) as a function of the internuclear 
separation r and of the energy Ec in the continuum is then obtained by fold- 
ing the corresponding transition density with the external field generated by 
the target. In addition to the Coulomb field, a Woods-Saxon nucleon-nucleus 
potential is used, with parameters of _R = vqA^^^, rg = 1.27 fm, a = 0.67 fm, 
V = (-51 + 33 (TV - Z)/A) MeV, and Vu = -0.44^. 

Selected cuts of the form factor F(r; Ec) are shown in Figure 1. The 
individual contribution arising from the Coulomb and nuclear interactions are 
shown separately. In Fig. 1 (a), we display the form factor as a function of 
r for a fixed value of the energy in the continuum {Ec being 0.9 MeV). Note 
the long tail of the nuclear contribution as a consequence of the large radial 
extension of the weakly-bound wave function, resulting in the predominance 
of the nuclear form factor up to the unusual distance of about 25 fm. The 
same reason gives rise to a deviation of the Coulomb part from the asymptotic 
behaviour 1/r^. Note also, due to the negative El effective charge of the 
neutron excitation, the constructive interference of the nuclear and Coulomb 
parts. In figs. 1 (b) and 1 (c), we show, instead, the energy dependence of the 
form factor at fixed value of r. While at large r the curves are peaked at very 
low energies, reflecting the corresponding B{E1), at smaller distances around 
the barrier, for values which are more relevant to the fusion process, the peaks 
of the distributions move to higher energies, in particular for the nuclear part. 
Hence, the sudden approximation employed in Refs. 7 and 8, where the energy 
of continuum states Ec was all set to be zero, is not justified. 

In order to perforin the coupled-channel calculation, the distribution of 
continuum states is discrctized in bins of energy, associating to each bin the 
form factor corresponding to its central energy. We have considered the contin- 
uum distribution up to 2 MeV, with a step of 200 keV. In this way, the calcu- 
lations are performed with 10 effective excited channels. We have checked the 
convergence concerning the maximum energy of the continuum states included 
in the calculation, and found that it is rather slow. We have, however, found 
that our main conclusions again do not qualitatively change provided that the 
maximum excitation energy of the continuum states is at least 2 MeV, as con- 
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Figure 1: Coupling form factor F{r;Ec) associated with the dipole transition in Be from 



the neutron bound state 2s 



1/2 



(Ei 



-500 keV) to the continuum state P3/2 with energy 



Ec- In (a) the Coulomb (the dash-dotted line), nuclear (the dashed) and total (the solid) 
form factors are shown as a function of r at the continuum energy Ec = 0.9 McV. In (b) 
and (c) the form factors are shown as a function of the energy Ec in the continuum for 
r = r[,£irrier=11.6 fm and r = 30 fm, respectively. 
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Figure 2: Fusion cross section for the reaction Be + Pb as a function of the bombarding 
energy in the center of mass frame. The thin soUd curve shows the results of the one- 
dimensional barrier penetration as a reference. The solid and the dashed lines are solutions 
of the coupled-channels equations for the complete fusion and the complete plus incomplete 
fusion, respectively. 



sidered in this work. The ion-ion potential is assumed to have a Woods-Saxon 
form with parameters Vq ~ —152.5 MeV, ro = l.l fm and a =0.63 fm, a set that 
leads to the same barrier height as the Akyiiz-Winther potential. 

In the previous theoretical studies of fusion of drip-line nuclei, LrB'EI the 
excitations to the "soft-dipole mode" of the projectile were also included in the 
calculations and the effects of the break-up were thus somewhat perturbed. In 
the present calculations, in order to isolate the genuine effect of the break- 
up process, we include only the continuum states in the coupling scheme, 
neglecting continuum-continuum coupling as well as other inelastic channels 
such as bound excited states in either reaction partner. For the same reason, we 
do not take into account static modifications on the ion-ion potential which may 
arise from the halo properties of the projectile. E3 Our calculations, therefore, 
only attempt to give qualitative indications. They will, however, still reveal 
interesting aspects of effects of couplings of the ground to continuum states on 
subbarrier fusion. 

Figure 2 shows the results of our calculations. The solid line represents 



the cross section of complete fusion, leading to ^^^Rn, while the dashed line 
denotes the sum of the complete and incomplete fusion cross sections. Also 
shown for comparison, by the thin solid line, is the cross section in the absence 
of the couplings to the continuum states. One can see that they enhance the 
fusion cross sections at energies below the barrier over the predictions of a 
one-dimensional barrier penetration model. Note that this is the case not only 
for the total (complete plus incomplete) fusion probability, but also for the 
complete fusion in the entrance channel. This finding supports the results of 
the original calculation performed in Ref. 9. As it has been emphasized there, 
accounting properly for the dynamical effects of the coupling in the classically 
forbidden region is essential to arrive at this conclusion. 

The situation is completely reversed at energies above the barrier. Fusion 
in the break-up channel becomes more important and dominates at the expense 
of the complete fusion. As a consequence, the cross sections for complete fusion 
are hindered when compared with the no-coupling case. 

3.2 The ^'^He + ^^^ U reactions 

We next consider the fusion reactions ''^'^He + 238tt, This is partly motivated 
by the recent measurement by the SACLAY group, u They compared the fusion 
cross sections of these two systems and concluded that the cross sections for the 
^He projectile were systematically larger than those for the '^He projectile, both 
below and above the Coulomb barrier. One might think that this conclusion is 
inconsistent with ours discussed in the previous subsection as well as with the 
experimental observation of the Canberra group u that the break-up hinders 
complete fusion cross sections at energies above the Coulomb barrier. However, 
the situation is aomewhat more complicated since the experimental data of the 
SACLAY group H likely include both the complete and incomplete fusion cross 
sections, and also since the halo structure of ^He may alter significantly the 
static fusion potential. We therefore decided to carry out coupled-channels 
calculations for these systems in order to clarify the role played by the break- 
up of the ^He projectile. 

To this end, we describe the structure of ^He using the di-neutron cluster 
model as in Ref. 23. The potential between the di-neutron and a particle is 
determined so as to reproduce the correct binding energy of ^He. The ground 
state is assumed to be the 2s state in this potential. In the coupled-channels 
calculations, we include continuum states with L = 1 and L = 2 spins. To 
construct these wave functions, we slightly adjust the depth of the potential for 
L = 2 so as to reproduce the correct resonance state at 1 MeV. We include the 
continuum states up to £^c = 5 MeV with energy step of 0.25 MeV. As for the 
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Figure 3; Same as fig. 2, but for the ®He + 238^ reaction. The meaning of each line is the 
same as in fig. 2 except for the dot-dashed fine, which denotes fusion cross sections for the 
4jjg _|_ 238u system. 



bare potential, wc use the Woods-Saxon potential with Vq = —174.05 MeV, 
ro = 1.01 fni, and a=0.67 fm for the "'He + ^38^ system. The potential for the 
®He -I- 238jj system is constructed with the folding procedure using the same 
di-neutron cluster model. For simplicity, we neglect the deformation effects 
of the 238^ target which would influence the both reactions in a similar way, 
although one definitely has to include them when one compares theoretical 
fusion cross sections with the experimental data. 

Figure 3 shows the theoretical fusion cross sections for the ^He + ^38^ 
reaction and their comparison with the ^He + ^sSjj system. The relation be- 
tween the complete fusion cross sections (the solid line) and the bare cross 
sections (the thin solid line) within the ^He + 238u system is qualitatively 
the same as that in the ^Be -I- 208pj^ system discussed in the previous sub- 
section; the coupling to the break-up channels enhances the complete fusion 
cross sections at energies below the barrier and hinders them at above the 
barrier energies. The cross sections for the *He projectile are denoted by the 
dot-dashed line. The figure shows that the complete fusion cross sections for 
the ^He -I- 238u reaction are even smaller than the cross sections for the ^He 
projectile at energies above the Coulomb barrier. Notice that the sum of the 
complete and incomplete fusion cross sections for the ^He projectile is always 
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larger than the cross sections for the ^He promctile, that remmisces the ex- 
perimental observation by the SACLAY group. lI It would thus be interesting, 
but perhaps experimentally challenging, to isolate the contribution from the 
complete fusion in the ^He + ^ssy gygtem. 

4 Summary 

We have performed exact coupled-channels calculations for weakly-bound sys- 
tems using realistic coupling form factors to discuss effects of break-up on 
subbarrier fusion reactions. As examples, we have considered the fusion of 
-^^Be with a ^°®Pb target as well as the '*'^He -I- ^38^ fusion reactions. We 
found that the coupling to break-up channels enhances cross sections for the 
complete fusion at energies below the Coulomb barrier, while it reduces them 
at energies above. 

Very recently, a complete fusion excitation function was measured for the 
^Be -I- ^"^Pb reaction at near-barrier energies by Dasgupta et al.a They showed 
that cross sections for complete fusion are considerably smaller at above-barrier 
energies compared with a theoretical calculation that reproduces the total fu- 
sion cross section. Also, the fusioii.cross sections for the ^'^He -I- ^38^ systems 
measured by the SACLAY group B seem to indicate that the break-up effects 
enhance fusion cross sections at energies below the barrier. These two recent 
experiments are in general agreement with our results. 
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